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We study laser cooling of two ions that are trapped in a harmonic potential and interact by 
Coulomb repulsion. Sideband cooling in the Lamb-Dicke regime is shown to work analogously to 
sideband cooling of a single ion. Outside the Lamb-Dicke regime, the incommensurable frequencies 
of the two vibrational modes result in a quasi-continuous energy spectrum that significantly alters 
the cooling dynamics. The cooling time decreases nonlinear ly with the linewidth of the cooling 
transition, and the effect of trapping states which may slow down the cooling is considerably reduced. 
We show that cooling to the ground state is possible also outside the Lamb-Dicke regime. We develop 
the model and use Quantum Monte Carlo calculations for specific examples. We show that a rate 
equation treatment is a good approximation in all cases. 



I. INTRODUCTION 



The emergence of schemes that utilize trapped ions or atoms for quantum information, and the 
interest in quantum statistics of ultra cold atoms, have provided renewed interest and applications 
for laser cooling techniques The present goal is to laser cool several atoms to a pure quantum 
state (to the motional ground state), and experimental [g] and theoretical efforts j| are made in 
this direction. The coolingof a large number of particles using lasers is a prerequisite for coherent 
control of atomic systems ji],^ . In quantum information, for example, laser cooling to the motional 
ground state is a fundamental step in the preparation of trapped ions for quantum logic . Coher- 
ent control and manipulation of information requires that each ion be individually addressable with 
a laser, and thus restricts the choice of the trap frequency, and consequently the regime in which 
cooling must work, to relatively shallow traps. 

Laser cooling of single ions in traps has been extensively studied Q and in particular sideband 
cooling has been demonstrated to be a successful technique for cooling single ions to the ground 
state of a harmonic trap Sufficient conditions for sideband cooling a two level system are: (i) 
the radiative linewidth 7 is smaller than the trap frequency v, such that motional sidebands, i.e. 
optical transitions that involve the creation or annihilation of a specific number of motional quanta, 
can be selectively excited; (ii) the Lamb-Dicke limit is fulfilled, i.e. the ion's motional excursion is 
much smaller than the laser wavelength. The first condition can be achieved through an adequate 
choice of atomic transition or a manipulation of the internal atomic structure while the Lamb- 
Dicke regime requires the trap frequency to be much larger than the recoil frequency of the optical 
transition. 

For more than one ion, as required in quantum logic schemes, individual addressing imposes small 
trap frequencies, whereas sideband cooling imposes high trap frequencies. Furthermore, the Coulomb 
interaction between the particles makes the problem much more complex, and it is not obvious 
whether the techniques developed for single ions can be transferred directly to this situation. Ex- 
perimentally, sideband cooling of two ions to the ground state has been achieved in a Paul trap that 
operates in the Lamb-Dicke limit . However in this experiment the Lamb-Dicke regime required 
such a high trap frequency that the distance between the ions does not allow their individual ad- 
dressing with a laser. For this purpose, and also for an extension beyond two ions, linear ion traps 
H are the most suitable systems, and for their physical parameters laser cooling to the ground state 
is a goal yet to be achieved. Laser cooling of two ions into the ground state is the problem that we 
address in this paper. 

Theoretical studies on cooling of single ions outside the Lamb-Dicke regime exist [TLlI 12|, while 



laser cooling of more than one ion has been analyzed only in the Lamb-Dicke regime [ 13|J14[ . In 
this paper, we investigate laser cooling, as developed for single ions, when it is applied to two or 
more ions. We focus our attention to sideband cooling, showing and discussing new physical effects 
which arise because of the presence of two interacting particles. Doppler cooling will be discussed 
in a future work. We show that cooling of two ions outside the Lamb-Dicke regime presents novel 
features with respect to single ion cooling, and we show how the preparation of the two ions in a 
pure quantum state is possible. The results will allow us to get some insight into the more general 
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problem of cooling a string of -/V ions. This is not only relevant for quantum logic with iV-ion strings 
but also for laser cooling of ion clusters in Paul and Penning traps ]6|Jl5| . 

The paper is organized as follows. In section II we introduce and discuss the model which we will 
use throughout the paper, and discuss some concepts developed in sideband cooling of one ion in 
relation to the presence of more than one ion. In section III we study and discuss sideband cooling 
of two ions inside and outside the Lamb-Dicke regime, and compare the two different behaviours. 
Finally in the conclusions we summarize the main results, and discuss the problem of cooling N > 2 
ions. 



II. MODEL 



We consider two ions of mass m and charge e placed in a one-dimensional harmonic potential 
of frequency v. We assume the ions to be strongly trapped in the other spatial dimensions so that 
their motion in those directions is frozen out. Their internal structure is described by a two level 
system with ground state |g), excited state |e) and resonance frequency wo- The ions interact with 
laser light at frequency ujl and wavevector k. For classical laser light and in the Rotating Wave 
Approximation the Hamiltonian of the system is: 

H = Hi + -f/mcc + V. 

Here Hi is the internal energy in the rotating frame 

Hi = -8 ^ c >j 

i=i,2 

where 8 — u)l — ujq is the detuning, j labels the ion (J = 1, 2) and we have taken % = 1. ff m0 c is the 
mechanical Hamiltonian, 
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with Xj, pj position and momentum of the jth ion (j — 1,2), and V describes the interaction 
between laser and atoms, 
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Here 0(x 3 ) is the Rabi frequency at the position Xj, crt, <jj are the raising and lowering dipole 
operators respectively defined on the jth ion (j = 1,2), and 8 is the angle between the laser 
wavevector and the trap axis. 

Using the center of mass (COM) and relative coordinates, the mechanical Hamiltonian in (^|) is 
composed of two separate terms: one for the COM motion which describes a particle of mass 
M = 2m interacting with a harmonic potential of frequency v\ the other for the relative motion which 
describes a particle of mass /i — m/2 interacting with a potential, which is the sum of a harmonic 
potential of frequency v and a Coulomb-type central potential. This potential may be approximated 
by a harmonic oscillator potential of frequency v r = obtained through the truncation at the 

second order of its Taylor expansion around the equilibrium distance Xq = (2e 2 /A-KeoMv) 1 ^ 3 between 
the ions |{l3|Jl4| . In appendix A we discuss this approximation and we show that it is valid in the 
regime that we are going to study. With this approximation the term (^[) becomes (apart for a 
constant) 

P 1 2 -i^2 P 1 2 2 

a mec = 1 — Mv X H 1 — nv r x , 

2M 2 2fi 2^ r ' 

where X — [xi + X2)/2, P = pi + pi are position and momentum of the COM respectively, and 
x = x\ — X2 — xo, p = (pi — P2)/2 are position and momentum of the relative motion. Thus the 
term apart for a constant, can be rewritten as 

Hmec = z^ojao + v r a\a r , 



where we have defined X = yJl/2Muo{a] } + ao)^ P = iyj Mvo/2{a\ l — ao), x = xq + ^l/2jiv r (a\.-\-a r ), 
p = iyj iiu r /2{a < r — a r ), with ao, aj the annihilation and creation operators for the COM mode 



respectively, and a r , a\. the corresponding ones for the relative motion (stretch) mode. We stress 
that in this new representation the mechanical problem of two ions interacting through Coulomb 
forces is reduced to the one of two harmonic oscillators, while the interaction of each ion with the 
radiation is now transformed into a nonlinear coupling between the harmonic oscillators. In general, 
N ions in a trap can be described by a set of N harmonic oscillators, coupled by laser light [jl7| . 
The master equation for the density matrix p of the two ion system is jl3 

-p = --[H, P } + Lp. 

Here L is the Liouvillian describing the incoherent evolution of the system: 

Lp = - ^2 l 2cr iPi cr ] ~ a \ a i? - P°\°i\, 

j = l,2 

where 7 is the decay rate out of the internal excited state |e), and pj describes the density matrix 
after a spontaneous emission for the jth ion: 

r 1 

p 3 =J duN(u)e lkux ipe- lkux i 
with N(u) being the dipole pattern for the decay. 

In this treatment we have neglected both dipole-dipole interaction between the ions and quantum 
statistical properties. This approximation is justified in the regime that we investigate, which is 
characteristic of experiments using linear ion traps for quantum information ||. In those traps the 
equilibrium distance between the ions is of the order of lOu, while the laser wavelength is typically 
in the visible region and the individual ionic wavepackets have spatial widths of the order of 10-100 
nm. From these considerations we can consider the two ions in a linear trap as two distinguishable 
particles in a harmonic potential which interact solely with Coulomb forces. On the basis of these 
considerations we will use Eq. (Q) for the numerical simulations presented below. 
For the following discussion it is instructive to look at the set of equations which one obtains from 
(Q) in the limit of low saturation Q <C 7, when the excited state |e) can be eliminated in second 
order perturbation theory (we provide a detailed derivation of the equations in Appendix B). In the 
basis of states |g, n), where n = (no,n r ) is a vector with COM vibrational number no and stretch 
mode vibrational number n r , we have the following set of equations: 



|p|g> m ) = -*( n - m ) ■ v is, n|p|g, m) 
\ " 



.n 2 



k,l 



(n|e ifccostel |k}<k|e- ifccosfla,1 |l) ^{l\e lkcosexi \k){k\e- ikcosexi \m) 



+ 



4 



(k-l)-u-<5-i7/2 > (k-l)-v -5 + 17/2 

V f d^( U )<n|c lfcM1 |k)<k|c- lfccostel |r)<s|c lfccostel |j)(j|e- l ^ 1 |m)<g,r|p|g, S ) 



+ 



1 
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where v = (v,l> r ) and where for simplicity we have assumed that only ion 1 is illuminated, i.e. 
Q(xi) = fl, Q.{x2) = 0. This implies that we can address the ions individually with a well focused 
laser beam. It corresponds, for example, to a situation where the two ions are two different isotopes, 
of which only one is resonant with light jl8| |. 

When treating laser cooling in a harmonic trap, an important dimensionless quantity is the Lamb- 
Dicke parameter rj, which for a single ion of mass m in a trap of frequency v, interacting with laser 
light of wavevector k is: 



1 / Wrec 

where ui TCC = k 2 /2m is the recoil frequency. The parameter n appears in the kick operator exp(ikx) 
in the term des cribing the exchange of momentum between radiation and atoms, which using the 
relation x = ^/l/2mw(a t + a) and the definition ([ll]) is rewritten as exp(ikx) — exp(iri(a^ + a)). 
The Lamb-Dicke regime corresponds to the condition \fnr\ 1, with n vibrational number; in 
other words, to the situation in which, during a spontaneous emission, a change in the vibrational 



number of the atomic state is unlikely due to energy conservation. In this regime the kick operator 
may be expanded in powers of 77, and with good approximation the expansion may be truncated at 
the first order Q|. Another important parameter, as known from cooling of single ions, is the ratio 
between the radiative linewidth 7 and the trap frequency v. in the so-called strong confinement 
regime 7/1/ < 1 the laser can selectively excite sidebands of the optical transition which involve 
a well-defined change of the vibrational number n. In this regime, together with the Lamb-Dicke 
regime, sideband cooling works efficiently: when the laser is red detuned with S — — v the system 
is cooled by approximately one phonon of energy v in each fluorescence cycle, finally reaching the 
vibrational ground state n = In contrast, in the weak confinement regime j/f > 1, transitions 
which involve different changes of the vibrational number n are excited simultaneously. This is the 
Doppler cooling regime, where the achievable minimum energy for a single ion is approximately 7/2 
for a detuning 5 = —7/2 Jl6| . 

Having introduced these basic concepts and methods of laser cooling of single ions in traps, we 
now turn back to the problem of two ions, to discuss how those techniques may be applied, and 
whether the same concepts are still valid. For N harmonic oscillator modes we can define a Lamb- 
Dicke parameter for each mode in an analogous way to (|ll|). For our case, N = 2, the Lamb-Dicke 
parameters 770 for the COM mode and rj r for the stretch mode are defined as: 



770 = k 



Vr = 7T 
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so that the kick operator for the jth ion (J = 1, 2) is written as: 

^ikxj _ e ivo( a l + a o) e (- 1 ) J ~ li Vr{al + a r ) 

In general for TV ions rjo = ry/x/TV jl?]]. In the following, when we refer to the Lamb-Dicke regime 
we will consider the situation where the conditions ^/no77o <C 1, ^/rwr\ r <C 1 are fulfilled. From @ 
we see that in Eq. (Q) the Lamb-Dicke parameters appear multiplied by the factor cos#. Therefore 
the Lamb-Dicke parameters for the coherent excitation, m cos 9, are always less than or equal to the 
Lamb-Dicke parameters defined in ([l]), which characterize the spontaneous emission. 
To discuss the importance of the ratio 7/V in the case of two ions, we first consider the bare spectrum 
of energies of our system with frequencies v and y?>v. In Fig. 1 we plot the density of states D(E) vs. 
the energy E, defined as the number of states in the interval of energy [E, E + SE] . From this figure 
we see that due to the incommensurate character of the frequencies the spectrum does not exhibit 
a well distinguished series of energy levels, rather it tends to a quasi-continuum. Therefore, the 
strong confinement requirement for sideband cooling needs to be reconsidered. The main question 
which we will address in the following is whether it is still possible to cool one mode to the ground 
state by means of sideband cooling. As we will show, the Lamb-Dicke parameter distinguishes two 
regimes which exhibit dramatic differences. 



III. SIDEBAND COOLING OF TWO IONS 



In the following we study sideband cooling of two ions, first in the Lamb-Dicke regime and then 
outside of this regime. We will show that in this latter case two-ion effects appear due to the dense 
spectrum of energy levels. In our calculations we first consider sideband cooling when laser light 
excites only one of the two ions directly. Thus in (Q) we take Fl(xi) = 57, and fl(x2) = 0. Afterwards 
we compare this case to the one in which both are driven by light, i.e. fi(xi) = £l(x2) = 0, showing 
that the only difference between the two cases is the cooling time, which in the latter case scales by 
a factor 1/2. In the following we assume the laser wavevector parallel to the trap axis, i.e. cos 6—1. 
This assumption facilitates the analysis and it is justified by the simple scaling just described. 
Furthermore it corresponds to the case in which the two ions are two different ionic isotopes, of 
which one is driven by light At the end of this section we will briefly discuss cooling of two 

identical ions when the wavevector is not parallel to the trap axis. 



A. Lamb— Dicke regime. 



In the Lamb-Dicke regime the Franck-Condon coefficients (n| exp(ifea;)|l) in the numerators of 
the right hand side terms of Eq. ( |Io| ) may be expanded in terms of the Lamb-Dicke parameters 770, 
rj r . The response of the system to laser light is governed by its spectrum of resonances 1(8), which 
is evaluated by summing all contributions to laser-excited transitions at frequency ujl, 

1(6) = ]T |<n|exp(fe)|l>| 2 ^(n) 

where P(n) is a normalized distribution of the states |n). In the Lamb-Dicke regime, we find 
that 1(5) exhibits two main pairs of sidebands around the optical frequency ljq: one at frequencies 
ljq ± vo corresponding to the transition no — > no ± 1; the other at frequencies uoo ± v r corresponding 
to n r — ■» n r ± 1 (see Fig. 2(a)). The strength of these sidebands relative to the carrier n — > n are 
proportional to 770 and to 77J:, respectively. All the other sidebands have strengths of higher orders 
in , rfc. This implies that by selecting one of these four sidebands by laser excitation we will 
induce the corresponding phononic transition; for example by choosing the sideband corresponding 
to no, n r — ► no — 1, n r we can cool the COM mode to its vibrational ground state, as for a single 
ion. This has been experimentally demonstrated by the NIST group at Boulder fllCf. In Fig. 2 we 
plot the results of a Quantum Monte Carlo wavefunction simulation (QMC) of Eq. (JtJ) for two ions 
in a trap with Lamb-Dicke parameter no = 0.1, radiative linewidth 7 = 0.2u, detuning 5 = — v and 
an initially flat distribution for the states with energy E < 15v. In Fig. 2(b) the average vibrational 
numbers of the COM mode (solid line) and of the stretch mode (dashed line) are plotted as a 
function of time in unit of fluorescence cycles tF — 27/fi 2 . The system behaves as if the two modes 
were decoupled, since only one mode is cooled while the other remains almost frozen. Nevertheless 
the stretch mode is cooled on a much longer time scale, as an effect of off-resonant excitation. In 
Fig. 2(c) the populations of the vibrational states of the two modes are plotted at time t = 600£f, 
showing the COM mode in the ground state and the nearly uncooled stretch mode. In this limit 
we can neglect the coupling of the population to the coherences in Eq. (|lo|), thus reducing the 
density matrix equation in the low saturation limit to rate equations. In fact the coherences in Eq. 
( px| ) have either an oscillation frequency which is much larger than the fluorescence rate 1/tF, or 
a coupling to the population which is of higher order in the Lamb-Dicke parameter expansion, or 
both of them. Therefore, they can be neglected in the equations of the populations, and we obtain 
the set of rate equations: 

d, Q 2 ^ l(n|e fel |k)| 2 . . . 

dt {nMn) = _7 T E [( k - n) .,-^ +7 74 <n|p|n) 
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4 L^J_ X K 1 [(k-r). W -<5] 2 + 774 V ' 

where we have omitted the label g of the states. The validity of this approximation is shown for 
the above case in Fig. 3(a), where the results of figure 2(b) are compared with those of a rate 
equation simulation according to Eq. (|l^). As a further proof that the two modes can be considered 
decoupled during the time in which the COM motion is cooled, in Fig. 3(b) we compare the time 
dependence of the average vibrational number of the COM mode with the one of a single trapped 
ion which is cooled under the same Lamb-Dicke parameter, radiative linewidth, trap frequency, 
Rabi frequency and initial distribution as the COM mode. We see that the two curves overlap 
appreciably, justifying the picture of sideband cooling of two ions in the Lamb-Dicke regime as if 
the modes were decoupled from one another. 



B. Outside the Lamb-Dicke regime 

To illustrate the physical features of the system outside the Lamb-Dicke regime, we plot in Fig. 
4(a) the spectrum of resonances 1(5) as defined in Eq. (0) for two ions in a harmonic trap with 
770 = 0.6. We see that the spectrum exhibits many sidebands whose density increases as the detuning 
increases. The main consequence is that we cannot select a given sideband by choosing the laser 
frequency, but rather excite a group of resonances that correspond to transitions to a set of quasi- 
degenerate states. The range of transitions that are excited increases with 7. In figure 4(b), (c) we 
consider sideband cooling for Lamb-Dicke parameter 7/0 = 0.6 and detuning 6 — ~2v W\\, where 



the other parameters are the same as in Fig. 2(b), (c). As one can see, the two modes are coupled 
and cooled together. Thus, as a first big difference with respect to single ion cooling, we see that 
here the energy is not taken away from one mode only, rather it is subtracted from the system as 
a whole. Another striking difference appears in the cooling time, which is significantly longer in 
comparison with the time necessary to cool one single ion outside the Lamb-Dicke regime |Tl| . This 
slowing down is partly due to the increase of the dimension of the phase space where the cooling 
takes place: the presence of two modes makes the problem analogous to cooling in a two dimensional 
trap, whose axis are coupled by the laser. The ions thus make a random walk in a larger phase 
space, and the cooling gets slower. However the cooling time is even considerably longer than one 
would expect taking the dimensionality into account. This can be explained by looking again at the 
spectrum of resonances in Fig. 4(a): despite the high density of resonances, the coupling between 
the states is still governed by the Franck-Condon coupling, i.e. by the terms in the numerator of 
Eq. (jil^) which outside the Lamb-Dicke limit oscillate with the vibrational numbers of the states. 
In the limit of linewidth 7 <C v, where a single sideband can be selected, we may encounter trapping 
states like in cooling of single ions ^|, i.e. states whose coupling to the resonantly excited state is 
very small since their motional wavefunction after the absorption of a laser photon happens to have 
a very small overlap with the motional wavefunction of the excited states. This effect limits the 
cooling efficiency, since the atoms may remain trapped in these states and not be cooled further, or 
much more slowly, towards the ground state. For two ions the probability of finding zeroes of the 
Franck Condon coupling is larger than for one ion, as the coupling to the excited state is constituted 
by two integrals, one for the COM and the other for the relative motion wavefunctions. Thus the 
probability of having trapping states is higher. To illustrate this phenomenon we plot in Fig. 5(a) 
the occupation of the states P n as a function of the COM and relative vibrational numbers no and 
n r , respectively, at a time t — 1000£f after sideband cooling of the COM with S = — 2v. Here 
7 = 0.02;/, and we are in the limit in which the single resonances are resolved. As a consequence the 
most likely coherent transitions are no — » no — 2, n r — > n r . The effect of the trapping states is visible 
in the tail of occupied states of P n , with n r = 6, 7. In Fig. 5(b), (c) we plot the modulus square of 
the Franck Condon coefficients for the relative motion corresponding to the coupling of the states 
n r = 6, 7 to the other motional states respectively: here it is clearly shown that for the transition 
n r = 6 -—> 6, n r — 7 — » 7, the coupling is reduced nearly to zero. As the linewidth 7 increases, the 
number of states to which a single state is coupled increases. Thus the number of channels through 
which the atom may be cooled is larger. As an effect the trapping states disappear. This is shown 
in Fig. 6, where the population P n is plotted for t = 600£f and 7 = 1/, and otherwise the same 
parameters as before. Here we see that the system is cooled homogeneously. The effect of varying 7 
is summarized in Fig. 7, where we compare the average COM vibrational number vs. time in unit 
tF for various values of 7. The results of Fig. 7 show clearly that as the linewidth increases the 
number of fluorescence cycles needed for cooling the system decreases dramatically. It is important 
to note that in this diagram the time is measured in units of fluorescence cycles for each 7, so that 
the absolute cooling time clearly reduces more strongly. We stress that this strong dependence on 
the linewidth is a two-ion effect. In contrast, in sideband cooling of single ions, the fluorescence 
time tF determines the cooling time scale for 7/1/ < 1, and the curves for different values of 7 vs. 
the time in unit of the respective tF do not show striking differences. 

The presence of trapping states and the coupling of each state to more than one state at almost the 
same transition frequency might lead to the formation of dark coherences between quasi-degenerate 
states, i.e. to superpositions of states which decouple from laser excitation because of quantum 
interference. However for the considered system those dark states do not play any significant role. 
We prove this numerically in Fig. 8, where we plot the comparison between a QMC and a rate 
equation simulation. We see that there are no striking differences between the two curves. We point 
out that outside the Lamb-Dicke regime a rate equations treatment is not justified in principle, 
since secular approximation arguments and Lamb-Dicke limit arguments cannot be applied. Here 
the rate equations are used to highlight the effect of neglecting the coherences in the dynamical 
evolution of the cooling, while these coherences are fully accounted for in the QMC treatment. In 
order to see why coherences do not play any significant role in the cooling dynamics, we look at 
the definition of a dark state. Let us consider a state |a) at t = defined for simplicity as linear 
superposition of two quasi-degenerate states \a) = ai[n) +a2C l ^|m), with ai, 02, 4> real coefficients 
and with n, m states almost degenerate in energy, so that |(n — m) • v\ = AE with AE < 7. In 
principle a dark state can be a linear superposition of any number of states. However as we will see 
from the arguments below our restriction to two states does not affect the generality of the result. 
The evolution \a(t)) in the Schroedinger picture, apart from a global phase factor, is written as 

\a(t)) =o x |n)+a 2 e < ' )iW |m> ) (16) 



6 



with 4>(i) — 4>o + AEt. The state is dark when the following condition is fulfilled: 

(\\e k(x+x,2) \a{t)) ~ 

for any state |1) belonging to the set of states {|1}} to which it is resonantly or almost resonantly 
coupled. If the condition (^) holds at t = 0, it will hold up to a time t such that AEt ~ tt/2. For 
the system we are dealing with we do not have exact degeneracy, thus we check whether the state 
| a) can remain dark for a time sufficiently long to affect the cooling dynamics appreciably. The 
smallest possible value of AE in the range of energies of our calculations is AE = 0.07f, and we 
find that 4> rotates by an angle 7r/2 in less than one fluorescence cycle for the values of 7 that we 
have considered. The dark coherences are then washed away during the evolution as an effect of 
the incommensurate frequencies between the two modes. This result together with the numerical 
results suggests that rate equations can be used in the study of cooling pfj . 

In order to highlight that the absence of dark coherences is a signature of the peculiar spectrum of 
the system we plot in Fig. 9 the cooling of one mode outside the Lamb-Dicke regime for the case 
of a discrete spectrum where we have exact degeneracy. More precisely we consider two harmonic 
oscillators with commensurate frequencies v and 2v, where all the other physical parameters are 
the same as before. In this case the different outcome between the QMC and the rate equations 
treatment is dramatic, giving evidence to the role of the coherences in the evolution. 



C. Light on both ions 

The calculations that we have shown refer to the case in which only one ion is illuminated. As 
we have seen, although light interacts with one ion it couples with both modes simultaneously, as 
shown in Eqs. (Q) and (|l3|). When both ions are excited by laser light, the system is described by 
a 4-level scheme, corresponding to the 4 internal states |ai,&2) with a,b — e,g, where we assume 
that when a photon is emitted, we detect from which ion the event has occurred, as a consequence 
of the spatial resolution of the ions. We expect that the effect on the cooling will be a doubling in 
the number of quantum jumps and hence of the cooling rate. This is shown in Fig. 10, where we 
compare the time dependence of the COM vibrational number for the cases in which only one ion is 
illuminated (dashed line) or both ions are illuminated with the same laser intensity (solid line with 
label 1). In the latter case cooling is visibly faster, and the time dependence scales with a factor of 
two with respect to the case with one ion illuminated, as we can see when we replot the solid line 1 
vs. t/2tf (solid line with index 2). The two-ions effects found above are clearly independent of the 
number of scattering points, with the only difference that the dark state condition is now written as 
(1| (fSHX+x/z) + e ik{X-x/2)^ ^ We point out again that e ff ec t s due to interference between 

the internal excitation paths have been neglected, as we consider the ions to be two distinguishable 
particles 

We would like to stress that in the above calculations we have considered the case of only one ion 
driven by radiation while the laser wavevector is parallel to the trap axis. However if one wants 
to cool two identical ions by shining light on one of them the laser beam must necessarily be at a 
certain angle 9 with respect to the ion string. Thus the Lamb-Dicke parameter characterizing the 
coherent excitation will be smaller than the Lamb-Dicke parameter for the spontaneous decay, and 
depending on the minimum angle 8 required, the coherent laser will excite with some selectivity 
one of the two modes. However for cooling purposes it is preferable to have the two Lamb-Dicke 
parameters values, corresponding to the spontaneous emission and the coherent excitation, as close 
as possible. 



IV. CONCLUSIONS 



In this paper we have studied the question of cooling two ions in a linear trap to the ground state 
by means of sideband cooling. We have studied sideband cooling in the Lamb-Dicke regime, and we 
have shown that in this limit the two harmonic oscillators can be considered decoupled when one 
of the two is cooled by means of sideband cooling. We have found that the cooling dynamics in the 
low intensity limit may be described by rate equations, and essentially that all the considerations 
developed for sideband cooling of single ions apply. This regime has been considered by Vogt et al. 
113] for studying the effects of dipole-dipole interaction in laser cooling of two ions and by Javanainen 
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in |l4| in his study on laser cooling of ion clusters. 

We have then investigated laser cooling outside the Lamb-Dicke regime, finding striking differences 
with the Lamb-Dicke limit. Here the energy spectrum may be considered a quasi-continuum, though 
the coupling between the states is still governed by the Franck-Condon coupling. A consequence is 
that the cooling efficiency depends strongly on the radiative linewidth. For very small ratio 7/1/ the 
effect of trapping states is appreciable, and manifests itself in a drastic increase of the cooling time, 
i.e. the number of fluorescence cycles needed to cool the system. For 7/v < f, but large enough, 
the effect of trapping states is washed away, the modes are cooled simultaneously, and the cooling 
time is considerably shorter and comparable to the time needed for cooling single ions outside of 
the Lamb-Dicke regime. These effects are all consequences of the density of states in the energy 
spectrum. A further property of the system is the absence of dark coherences, as a consequence of 
the incommensurate frequencies of the harmonic oscillators, i.e. of the absence of perfect degeneracy. 
This implies that in the low intensity limit rate equations provide still a good description of the 
cooling dynamics. Finally we have compared cooling when light is shone on one ion only or on both 
ions, finding a simple difference of a factor two in the rate of cooling. 

On the basis of the obtained results we would like to comment on cooling of a string of TV > 2 
ions. TV ions in a harmonic trap may be described by TV harmonic oscillators. The mode frequencies 
1/0, v\, vn-i are all incommensurate, and the number of states in the interval of energy [E, E+8E] 
is D(E) = E /voi/i...i/ff-iSE. Outside the Lamb-Dicke regime the spectrum of resonances is 
then even more dense than in the two-ions case, and the probability of trapping states will be larger. 
However a suitable increase of the linewidth will cancel their effects, since each state will see an even 
higher number of states, thus of possible transitions, than for two ions. We expect that this last 
property will also play a role against the creation of dark states. In fact, although on one hand the 
large density of states may favour the appearance of stable dark coherences between "accidentally" 
degenerate states, on the other hand for values of 7 large enough the coherent effects will wash 
out because of the coupling to a "continuum" of states. From these considerations we expect laser 
cooling to the ground state to be still possible for TV > 2 ions. Furthermore, one can cool a given 
set of modes to the ground state through the choice of the laser detuning outside the Lamb-Dicke 
regime. For example taking a detuning S — —v^ the modes with frequency Vj > may be cooled 
to their vibrational ground states. It should be noticed that as the number of ions TV increases, the 
Lamb-Dicke parameter of each mode decreases approximately as 1/VTV |17| , allowing to reach the 
Lamb-Dicke regime also when this condition is not fulfilled for single ions. In this limit the modes 
may be considered decoupled and sideband cooling is particularly efficient. 

As a last consideration, we note that the behaviour of a number of ions TV > 3 cooled by light 
depends on which ions of the string are driven. In fact each position of the string couples with 
the different modes with amplitudes that depend on the position itself ]T?j ]. Only if all the ions 
are illuminated we may consider all the modes as coupled and cooled simultaneously. But this 
"coupling" changes as we select and drive only certain ions of the string, fn this case a certain 
amount of modes may be cooled, while the others will remain hot or get cooled on a longer time 
scale. In this respect the system can be considered as having a reduced dimensionality, and the time 
of the cooling will be accordingly shorter with respect to the case in which the laser couples to all 
modes. 
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APPENDIX A: HARMONIC APPROXIMATION 

We consider the term (^) and rewrite it in COM and relative motion canonical variables: 



P 2 , 1„ 2^2 , P 2 , 1 



2M 2 2p 2 r 4ne \x\ 

H(X,P) + H(x,p). (AI) 
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with /i = m/2 reduced mass, M = 2m total mass. The mechanical problem is separable into center 
of mass motion and relative motion, where H(X, P) describes the harmonic motion of a particle of 
mass M interacting with a harmonic oscillator of frequency v, and H{x,p) the motion of a particle 
of mass (i interacting with a potential V(x) — fiu 2 x 2 /2 + e 2 /4ty€o\x\, i.e. a harmonic potential of 
frequency v and a central repulsive Coulomb-type potential. 

We focus our attention on the potential V(x), restricting its domain on the semiaxis x > 0. The 
equilibrium point is found to be xq = (e 2 /47reo/w 2 ) 1 3 . Expanding V(x) around xq we find: 



v{x) = \ {-^r) + \^ 2{x ~ xo)2 + T, { ~ ir i^xY i{x ~ xoT 



n— 3 



' ' -/if 2 J + \liv 2 (x — xo) 2 + A{x), 



2 \47reo / 2 

where v r = y/3u and A(x) sum over the higher order terms, which we call the anharmonic terms. 
We quantize the oscillation around the equilibrium position xo; 



xo + 



where a r , at are annihilation and creation operators respectively. From perturbation theory we may 
consider V(x) harmonic when the following conditions are fulfilled: 

(j\A(x)\j) < h »r 

\{j\A(x)\j±i)\<:nvr 

with \j) eigenstate of the harmonic oscillator of frequency v r with eigenvalue jhu r . The first con- 
dition means that the energy shift due to the anharmonic term is much smaller than the spectrum 
separation, whereas the second condition means that the coupling between the states is a small 
perturbation, and we will show that it may be neglected. The coupling between the state j and the 
state j + k is not taken here into account for simplicity, but it may be shown that it is m uch smaller 



than khv r in a similar way to the one we discuss below. Let us rewrite the relations in ( A5 ) as 



(j\A{x)\j) = 

4-KfnTn 



4-ireoXo 

2 00 

^^■ + 1 > = -4^£ C ^ +1 

m— 1 

with 

\ 2m 

1 

X^lpuVr 

where Cj,2m+i is analogously defined. From the following relation 



0.2m = ( \ j „ 2 ) 01 (at + a r ) 



i ,.,/ + \ 2m , -\ (2m)! (j + m) 



,2m 



< ¥ 



J'- J- 

i , , / t \2m+i. , (2m +1)! (j+rn)\ , 9™+i(? + m)! / 

y/j + m + 1 < (j (at + a r ) j ± 1 < \ U '- + m + 1 < 2 2m+1 u ' + ; 

v v ' m!(m + 1)! r- J- 



wc find that 




From ( |A9[ ) we see that the series does not converge, as the term of the expansion depends on the term 
m. However for a certain interval corresponding to j, m < Mo such that ip — ^Jl/ x 2 2[iv r ^/ Mo <C 1, 
each term is bounded by the corresponding term of a geometrical series with factor ^ < 1. For 
typical values of a linear ion trap |j| ip ~ y^Mo/100, so that Mo can assume very large values, 
Mo < 10 4 . 



The divergence at orders m > Mo is a signature of the divergence of the Coulomb potential at 
x — 0: such divergence requires the wavefunctions to be zero in x < 0, whereas the harmonic 
oscillator wavefunctions are different from zero on all the space, although their occupation on the 
semiaxis x < is very small for ijj <C 1. The divergence at x — is a mathematical consequence, 
and does not correspond to the physical situation, in which the ions move in a three dimensional 
space, although the confinement in the other two dimensions is relatively tight. In principle one 
could build up a potential that does not diverge in a specific point of the x axis, whose behaviour 
for x > tends to the Coulomb one. In this limit the series should converge for low number states 
j. However such a detailed study is beyond the scope of this paper pj^ |. We restrict to the case in 
which tp <C 1, showing that in the chosen range the harmonic approximation is quite good. 
From this consideration it is then sufficient to compare the third order term with the harmonic 
potential in order to show that the approximation is sensible. The relation to be satified is then 

%Vr > a 6 4 (\hr~l A 

47re a;o \ V 2 ^ / 

with j < Mo, which poses a further condition on j, Mq. Manipulating the expression we find 

_ / Hu r \ ( XQ \ 

J max ~ y e 2/ 4ne(jXo J y v /n/2flV r J 

which substantially agrees with the qualitative estimate in ^2|. For j j ma x the potential can be 
considered harmonic. For linear traps |(| j max ~ 120, and the harmonic approximation is valid for 
the region of energy we are considering. As an example the case j — 100 corresponds to a correction 
of the order of I0~ 3 hv r . 



APPENDIX B: ADIABATIC ELIMINATION OF THE EXCITED STATE 

We rewrite (Q) in the following way 

d i + 

&t p = ~n [HeSp ~ 9 ^ + Jp ' 

where H c g is the effective Hamiltonian 

H e s — Hi+ Hmec — crjcTi 



2 

i=l,2 



and Jp the jump operator 



We introduce the Liouvillians: 



Jp = 7 2J cnpaj 



L lP = -^[V,p] 
L2P = Jp 



so that (Bl) can be rewritten as: 



^p = [Lo + Li + L2] p 

In the limit Q <C 7 we can eliminate the excited state in second order perturbation theory. Calling P 
the projector onto the internal ground state |g), and using a standard derivation based on projectors 
H], we obtain the following equation for the ground state of the system: 

—Pp{t) = PL Pp(t)+ / d T PL 1 (l-P)exp(L T)(l-P)LiP P {t-T) 







dt 

+ j dn f dr 2 PL 2 (l-P)exp(L T 1 )(l-P)L 1 exp(Lo(r 2 -r 1 ))Pp(t-r 2 ) 

JO J T\ 
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where P is a projector so defined on a density operator X PX = g)(g|A|g)(g|. Markov approxima- 
tion can be applied in the limit in which we may consider the coupling to the excited state to evolve 
at a higher rate respect to the time scale which characterizes the ground state evolution. This is 
true once we have moved to the interaction picture with respect to the trap frequency. We define 

V](t) = e iBmtct Pp(t)e~ iHmm:t 



and in the interaction picture (B7) has the form 



vi (t) =e l 



drPLi exp {Lot) Lie 



-i-ffmoc(t-T) 



vi (t — r)e 



iffmec(*-T) 



+ c ■ 



dri / (IT2PL2 exp (ion) L\ exp (Lo(j 2 — n)) CLie 



-iHmcc(t-T 2 ) 



Vi[t - r 2 )e 



Now we can neglect the change of vi during the time r on which the excited state evolves. Going 
back to the original reference frame we have now the equation in Markov approximation: 



Pp(t) 



drPLi exp (L r) L l e lH ° T P p{t)e 



-m T 



+ dn dr 2 PL 2 eiqy{L T 1 )L 1 ^p{L (T 2 ~r 1 ))CL 1 e tHo ^Pp{t)e- lHoT2 
Jo J T1 



We substitute now in (B1C) the explicit form of the operators. After some algebra and application 
of the commutation rules we obtain the following equation (where we have neglected the interference 
terms between the two ions): 



d p(t)=PL 0P (t)-j2^-p 



dt 



+ e 



i=l,2 
-(-iA+ T /2)r 



i=l,2 



4 



4 

dri / dr2 / duN (u)e kuXi 



\ e 'e e e 



-iH m rcT ikx; iH mr .cT — ikx: 
moc g ,g moc g , 

r-1 



[e 



TL—Tz) //2(7"i+T 2 ) -iHmecTl -ifaj iilmecTl , 



<• '*•"" '«• •e"-""»" l p 1 (i)e """ """■■'e'" *e" 



met. ^„ z„ met. z 



+iA(-ri-T 2 ) — f/2(r 1 +T 2 ) -iH Blec T 2 -ikxi iH IaBC T2 



-ikux.i 



Projecting (Bll) on the basis of states {|n)}, we obtain Eq. (O) 



(B8) 



(B9) 
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FIG. 1. Density of states D(E) plotted as a function of energy E in units of v. The grid is SE — v /3. 



FIG. 2. (a) Spectrum of resonances for no = 0.1 for a thermal distribution with average energy per mode nu — 7.5v, plotted 
on a grid of width v/10. (b) Plot of (no) (solid line) and (n r )(dashed line) vs. time in units tF = 2^/Q 2 , for no = 0.1, 7 = 0.2v, 
Q = 0.034^, S = —u, and atoms initially in a flat distribution on the states with energy E < 15v. (c) Population of COM mode 
P n<) (onset) and of stretch mode P„ r (inset) vs. the respective vibrational state number at t — 600£f- 



FIG. 3. (a) Comparison between rate equation (solid line) and QMC calculation (dashed line). Same parameters as in Fig. 
2(b). (b) Comparison between the time dependence of the COM average vibrational number as in Fig. 2(b) (solid line) and 
the average vibrational number for the case in which a single ion is cooled (dashed line). For the single ion the mass has been 
rescaled so that = no, =!/,) = 0.2^, Q, = 0.034«/, with an initially flat distribution for the first 15 states. 



FIG. 4. (a) Spectrum of resonances for no = 0.6 of ions in a thermal distribution with average energy per mode nu = 7.hv. 
(b) Plot of (no) (solid line) and (n r )(dashed line) vs. time in unit tF = 27/fi 2 , for no = 0.6, 7 = 0.2u, Q. — 0.034^, 5 = —2v, 
and atoms initially in a flat distribution for states with energy E < 15v. (c) Population of COM (onset) and of relative motion 
(inset) vs. the respective vibrational number state at t = 6004f. 



FIG. 5. (a) Population P n as a function of no and n r at a time t — 1000£f and for 7 = 0.02^, Q = 0.177 = 0.0034^, 
tF = 27/fi 2 = 3460/za All the other physical parameters are the same as in Fig. 4(b), (c). (b),(c) Modulus square of the Franck 

Condon coefficients for the relative motion Fi } „ r = |(Z|e l ' 7 '"' £Ir+0 '"' |n r )| 2 with I = 6,7. 
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FIG. 6. Population P n as a function of no and n r at a time t = 600iF and for 7 = v , Q = 0.177 = 0.17u, t F = 2-y/Q, 2 = 69/za 
All the other physical parameters are the same as in Fig. 4(b), (c). 



FIG. 7. Time dependence of the average vibrational number of the COM mode for 7 = 0.02u (dashed line), 7 = 0.2v (solid 
line), 7 = OAv (dashed-dotted line) and 7 = v (dotted line), keeping constant the ratio fi/7 = 0.17. The time is in unit 
£f(7) = 27/fi 2 w 7O/7. All the other parameters are the same as in Fig. 4(b), (c). 



FIG. 8. Comparison between rate equation (solid line) and QMC calculation (dashed line). Same physical parameters as in 
Fig. 4(b). The onset refers to the COM and the inset to the relative motion vibrational number. 



FIG. 9. Plot of the average vibrational number vs. time for the harmonic oscillator of frequency v coupled to a second one 
with frequency 2v. Comparison between rate equation (solid line) and QMC calculation (dashed line). Lamb-Dicke parameter 
for the mode v \i, r\ v — 0.6, 7 = 0.2i/, fi = 0.034^, S = —2v, and atoms initially flat distributed on the states with energy 
E < Vov. 



FIG. 10. Plot of the average vibrational number of the COM mode vs. the time in unit tF for the case in which both ions 
are illuminated (solid line, index 1) and only one ion is illuminated (dashed line). The solid line with index 2 corresponds to 
line 1 rescaled, where the time has been divided by a factor 2. The other parameters are the same as in Fig. 4(b), (c). 
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